HEAVY QUARK PRODUCTION IN 77 COLLISIONS 



JIRf CHYLA 

Center for Particle Physics, Institute of Physics, Academy of Science of the Czech 
Republic, Na Slovance 2, 18221 Prague 8, Czech Republic 
E-mail: chyla@fzu.cz 

Current theoretical framework for the calculation of heavy quark production in 77 
collisions is reviewed. The importance of including direct photon contributions up 
to the order o?o? a and of the proper choice of renormalization and factorization 
scales in the evaluation of (7(77 — > QQ) is emphasized. 



1 Introduction 

Heavy quark production in 77 collisions has recently received increased theo- 
retical attention Ercra motivated in part by new experimental data on cc and 
bb production from LEP2 (see □ for references). Although the data on bb 
production have sizable errors and the theoretical predictions suffer from un- 
certainties, the excess by a factor of about three of the data over theoretical 
predictions represents serious problem for perturbative QCD. Interestingly, 
similar excess has been observed also in 7*p and pp collisions. 

In such a situation it appears timely to reanalyze the theoreticalpirame- 
work used for analyses of heavy quark production in 77 collisions pay- 
ing particular attention to the renormalization and factorization scale depen- 
dence, as these represent the main source of the theoretical uncertainty. 

2 Basic facts 

The factorization scale dependence of parton distribution functions (PDF) of 
the photon is determined by the system of coupled inhomogeneous evolution 
equations for quark singlet and nonsinglet and gluon distribution functions 



6S(x,M) 
dlnM 2 

dG(x,M) 
dlnM 2 
dq NS (z,M) 

dlnM 2 



= 5^k q + P qq ® E + P qG <g) G, (1) 
= k G + P Gq ® E + P GG ® G, (2) 
= S NS k q + Pns ® <7ns, (3) 



where, for n/ massless quark flavors, E(ir,M) = ^2™=i(qi(x, M) + q^x^M)), 
Ss = 6n/(e 2 ) and qjqs(x,M) and 5ns can be defined in two different ways. 
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The splitting functions admit perturbative expansion 



k q (x,M) 



2n 



*^) + ^^(*)+(^) a *«(*)+--- 



k G {x,M) = — 



Pij{x,M) 



a s {M ) . 
q.(M) 

2tt 



where k^ (x) 



(x 2 



(l-xf ) and 4 0) (x) are unique, whereas all higher 

(i) 

H )j > 1 depend on the choice of the factoriza- 
tion scheme. The equations can be recast into evolution equations 
for qi(x, M), q^x, M) and G(x,M) with inhomogeneous splitting functions 



a s (M) 

a,(M) 
2tt 
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3e 2 kq°\ The couplant a s depends on the renormalization scale fi and 



satisfies the equation 
da a (/x) _ 



/?Mm)) 



dln/i 2 

where for rif massless quarks /3q 



A) 2/ \ 



ft 
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11 - 2W3 and /?i = 102 - 38ra//3. 
General solution of the evolution equations (^Q) can be written as the sum 
of a particular solution of the full inhomogeneous equations and a general 
solution, called hadron-like (HAD), of the corresponding homogeneous ones. 
A subset of the former resulting from the resummation of contributions of 
diagrams describing multiple parton emissions off the primary QED vertex 
7 —>■ qq and vanishing at M = Mq, are called point-like (PL) solutions. Due 
to the arbitrariness in the choice of Mq the separation 

D(x, M) = D PL {x, M, M ) + L> HAD (x, M, M ). (8) 

is, however, ambiguous. Because of a different nature of the UV renormaliza- 
tion of the couplant a e (ri and the mass factorization involved in the definition 
of dressed PDF I will keep the renormalization and factorization scales /i and 
M as independent free parameters. 

In UnB the "next-to-leading order" QCD approximation to 17(77 ~~ * QQ) 
is defined by taking into account the first two terms in the expansion of direct, 
as well as single and double resolved photon contributions 

- «7<£ (M, fi)al (ja) + (M, rial (/x) + ■ 
a£ 2 ) (M, ri)a 2 s {ri + <4 3) (M, ri<{ri + ■ ' 



(0) 



a^(M)a s (ri 



',(9) 

•,(io) 



a^(M)at(ri + o-£'(M,rial(ri + • - (H) 



,(3) 
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Three light quarks were considered as intrinsic for cc and four in the case of 
bb production. Note that cr^ = cr c( ) (s / tUq) , where Co = 127ra 2 eg/s and 
c-°'{x) is a function of s/voq only, comes from pure QED. 




approach this is justified by claiming that PDF of the photon behave as a/a s 
and, consequently, all three expansions in fact start and end at the same 
powers (a s )° = 1 and a s , respectively. However, as argued inO, the logarithm 
In M 2 characterizing the large M behaviour of PDF of the photon comes from 
integration over the transverse degree of freedom of the purely QED vertex 
7 — > qq and cannot therefore be interpreted as l/a s (M). If QCD is switched 
off by sending, for fixed M and Mq, A — > 0, the quark and gluon distribution 
functions of the photon approach the QED expressions 

qi (x,M) - 9 ? ED ee £ 3ei 2 fc(°)(x)lngJ, G(x,M) - G^ ED = 0. (12) 

This is manifestly true for the point-like parts of quark and gluon distribution 
functions, whereas the vanishing of the hadron-like parts g HAD and G HAD as 
A — > follows from the fact that vector mesons (or in general any hadronic 
structure) can develop from the photon only as a consequence of strong in- 
teractions between the qq pair coming from the primary, purely QED, vertex 
7 — > qq. As QCD coupling vanishes, so must these hadronic parts. Let me 
emphasize that there is no obstacle to performing this limit, as by decreasing 
A we get ever closer to the asymptotic freedom point a s — and thus our 
perturbation expansions are progressively better behaved. In the limit of van- 
ishing colour coupling we thus get the purely QED result cr^ . Had the PDF 
of the photon really behaved as a/a s , we would get finite contributions from 
the lowest order single and double resolved photon contributions (|l(|[ll]) even 
in the limit of switchiufr.QCD off, which is clearly untenable. 

All calculations mm have been performed with fixed quark masses, i.u. 
define ttiq in the on-shell scheme. In this convention a^ t = <j$c^ (s/mg), 
where c^(x) is again a function of s/voq only. 

3 Direct photon contribution 

For proper treatment of Odin the total cross section for e + e _ annihilations 
into hadrons provides a suitable guidance. For n/ massless quarks we have 

f J had(Q) - 4ad + « S (M)4ad(Q) + «2 OO^Sd (Q / A + " " = 4ad(l + KQ)), 

(13) 
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Figure 1. Examples of diagrams describing the direct photon contribution to 17(77 ~* QQ) 
up to the order o?o? B . The solid (dashed) lines denote light (heavy quarks). 



where = {Aira 2 / Q 2 ) J2f e / comes, similarly as in (|9|), from pure 

QED, whereas genuine QCD effects are contained in the quantity 



(0) 



1 1 " s ^) in 1 \ 
1 H r^Q/fi) 



(14) 



In the case of the quantity (|T^) nobody calls the lowest order term <Jj^ d (Q) 

the "LO" and the next one, i.e. a^^ A {Q)a s {p)/'K, the "NLO" QCD approxi- 
mations, but these terms are reserved for genuine QCD effects described by 
r(Q)\ To work in a well-defined renormalization scheme (RS) of a s requires 
including in ([li]) at least the first two (nonzero) powers of a s (fj,) because the 
dependence of the coefficients rk(/J,/Q) on the RS starts with r\. The explicit 
/i-dependence of r\{Qj \x) cancels to the order a 2 s the implicit renormalization 
scale dependence of the leading order term a s (/i)/7r in (|l4| ) and thus guar- 
antees that the derivative of the sum of first two terms in (|TJ) with respect 
to In \x behaves as a^.. Because a s (/i) is a monotonous function of /i span- 
ning the whole interval (0,oo), the inclusion of first two consecutive nonzero 
powers of a s is also a prerequisite for the applicability of any of the scale 
fixing methods (seeS for detailed analysis of renormalization scale dependence 
of ([li])). For purely perturbative quantities the association of the term "NLO 
QCD approximation" with a well-defined renormalization scheme is a gener- 
ally accepted convention, worth retaining for physical quantities in any hard 
scattering process, like the direct photon contribution o^ir in ([§)■ 

Unfortunately, in the QED contribution a^?) is considered as the LO 
and the sum of the first two terms i n (p| ), which for fixed uiq has exactly the 
same form as the first two term in (|f3Q, as the NLO approximation. Conse- 
quently, this approximation cannot be associated to a well-defined renormal- 
ization scheme of a s and therefore does not deserve the label "NLO" even if 
the NLO expression for a s (fj,) is used therein. For QCD analysis of Udir in a 
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well-defined renormalization scheme of a s ([i) the incorporation of the third 
term in (0), proportional to a 2 a 2 , is indispensable. 

At the order a 2 a 2 , diagrams with light quarks appear as well and we can 
thus distinguish three classes of direct photon contributions differing by the 
overall heavy quark charge factor CF: 

Class A: CF = eg. Comes from diagrams, like those in Fig. [l]a-g, in which 
both photons couple to heavy QQ pairs. Despite the presence of mass 
singularities in contributions of individual diagrams coming from gluons 
and light quarks in the final state and from loops, at each order of a s the 
sum of all contributions of this class to o^ir is finite. 

Class B: CF = eg. Comes from diagrams, like that in Fig. |l|h, in which one 
of the photons couples to a heavy QQ and the other to a light qq pair. 
For massless light quarks this diagram has initial state mass singularity, 
which is removed by introducing the concept of the (light) quark and 
gluon distribution functions of the photon. 

Class C: CF = 1. Comes from diagrams in which both photons couple to 
light qq pairs, as those in Fig. 01. In this case the analogous subtraction 
procedure relates it to the single resolved contribution of the diagram in 
Fig. ||f and double resolved contribution of the diagram in Fig. ||h. 

As the diagrams in Fig. [j]e and [l]L give the same final state qqQQ, we have 
to consider their interference term as well, but it turns out that it does not 
contribute to the total cross section (7(77 — > QQ). Because of different charge 
factors, the classes A, B and C do not mix under renormalization of a s and 
factorization of mass singularities. At the order a 2 a 2 all three classes of 
contributions are needed for theoretical consistency. As argued above, class 
A is needed if the calculation of Odir is to be performed in a well-defined RS. 



4 Resolved photon contribution 



The classes B and C of direct photon contributions are needed to render the 
sum of direct and resolved photon contributions factorization scale invariant. 
To see this, let us write the sum of first two terms in ( |io[|ll| ) explicitly in 
terms of PDF and parton level cross sections 



T (12) 



(M, n) = 2a s (n) dxG(x, M) 



t (1) 



(X) 



(2) 



ia 2 M dx M)o$ (x, M) 
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Figure 2. Examples of diagrams describing resolved and related direct photon contributions. 



2a 3 a (ji) J jdxdy^x, M)G(y, M)ag (xy, M) + 

2a 2 s ( l i)JJdxdyJ2<li(x,M)q i (V,M) [<r^{xy) + a.(ji)a^{xy,M,ii)] + 
a 2 s ( f i)ffdxdyG(x 1 M)G(y,M) \a ( G ] G (xy) + a s {fi)af G {xy, M, n)] , (15) 



where % runs over n f quark flavors and the factors of two and four reflect 
the identity of beam particles and equality of contributions from quarks and 
antiquarks. Recalling the general form of the derivative dcr res /dlnM 2 



dov 



dlnM 2 



da; W (x, M) + dx 



dxdy 



£ q^(x, M)W qi (x, M) + G(x, M)W G {x, M) 



G(x, M)G(y, M)W GG (xy, M)+^ qi (x, Mfofo, M)W q - q (xy, M) 



+E(x, M)G(y, M)W qG (xy, M)\ , 
using and denoting a s = a s (/i), f = d//dlnM 2 we find 

Wo = ^{ k k ° {x)(J ^ {x) + 6fc * 0) {x) E 4°% & M ) } + • • • 

W qi = ^(W 2 )(x) + jdz\p^ Q (z)a^(xz) +Ze?akf\z)a ( ji\xz) 



W G = ^2^(x)+fdzP^(z)a^(xz)^ + 
W GG = ^ ^ G (x) + JdzpW(z)v, 



gU xz ) }+■ 



(16) 

(17) 

(18) 
(19) 
(20) 
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W qq = ^ l2na^(x) +2jdzP^(z)a^(xz)^ + ■■■ (21) 
W qG = ^ {2™^) + Jdz[p^(z)a^(xz) + Pg}{z)*™(xzj\}+- ' '(22) 



All integrals in ( |18||22j ) go formally from to 1, but threshold behaviour of 
cross sections <?ij(xz) restricts the region to xz > Aitiq/S. The factorization 
scale invariance of ( |l5| ) requires that the variation of ( |l5|) with respect to 
InM 2 is of higher order in a s than the approximation (M) itself. There is 
no question that direct photon contributions of classes B and C are needed 
to make the sum direct and resolved photon contributions factorization scale 
independent. The difference between the conventional and my approaches to 
QCD analysis of (7(77 — * QQ) concerns the question when this happens, i.e 
what is the order of the approximation ( |l5| ) and, consequently, which terms 
on the r.h.s. of (|l^ ) must vanish. 

In the conventional approach both q(M) and G(M) are claimed to be of 
order a/a s , and the approximation (|l5| ) thus of the order a 2 a s . This implies 
that only terms up to this order must vanish in (|l6|). As the expressions 
in round brackets of (|l8| ) and ([Ts| ) do, indeed, vanish, the expression (15) is 



claimed to be complete NLO approximation to cr ros . The fact that the expres- 
sion on the r.h.s. of ( |l7| ) docs not vanish is of no concern in this approach as it 
is manifestly of the order aa 2 and thus supposedly of higher order than ( |l5| ) 
itself. This line of argumentation is, however, untenable on physical grounds 
as it requires that not only the quark but even the gluon distribution func- 
tion G{M) should behave as a/a s , despite the fact that gluons appear in the 
photon only due to their radiation off the quarks and their contribution to 
physical observables must therefore vanish if this radiation is switched off! 

If we discard the untenable claim that q, G cx a/a s and realize the obvious, 
namely that q(M) cx a, we are lead to the conclusion that Wo is of the same 
order a 2 a 2 , as qiW qi and other integrands in (jig ) and must also vanish for 
theoretical consistency of the approximation ([151). This, in turn, necessitates 
the inclusion of class B direct photon contributions of the order aa 2 , like those 
in Fig. ^b,g, which provide the A/-dependcnt terms the derivative of which 
cancel Wq in the expression for d(a- ( a I (M) + fT res (Af))/dlnM 2 . Note that W qi 
in ( |l8| ) receives contributions from the derivatives of both single and double 
resolved photon diagrams, proportional to <t 7 g and a qq , respectively. This fact 
reflects the mixing of single and double resolved photon contributions, which 
starts just at the order a 2 a 2 . Note also, that for theoretical consistency of 
the sum Udu- + c ros up to this order, only the lowest order double resolved 
photon contribution needs to be included. 
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5 Phenomenological implications 



Although the current calculations of er(77 — > QQ) do not represent complete 
NLO QCD approximation, one may ask how important are the missing terms, 
discussed above, numerically. To help answer this question it would useful if 
the data could be expressed in terms of energy dependence of (7(77 — > QQ). 
The point is that for bb production direct photon contribution dominates 
below about y/s ~ 25 GeV, whereas for 200 > y/s > 50 GeV single resolve 
photon contribution makes up most of (7(77 — > bb). If the disagreement comes 
from the former region, we would be in real trouble. As, however, the LEP 
data integrate over the whole WW spectrum, to which the direct and single 
resolved photon channels contribute roughly equally, there is no way how to 
determine wherefrom comes the disagreement withlhe data. 

As argued above, the approximation used in naa for the evaluation of the 
direct contribution is merely of the LO nature. On the other hand, as the 
lowest order QCD correction crQa s makes up a small fraction of the purely 
QED contribution , the eventual incorporation of the (so far unknown) 
direct photon term crS^a^, coming from class A direct photon diagrams, will 
unlikely enhance <7di r by more than further few tens of percent. n 

The situation is different for the resolved photon contributions. In u I will 
discuss numerical aspects of the factorization and renormalization scale depen- 
dence of the approximation dl5| ) and show that proper choice of these scales 
is essential for understanding the data on bb production and may explain part 
of the observed discrepancy between the data and theoretical calculations. 
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